Structural Signatures of Vitrification in Hard Core Fluids 
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Abstract 

Computer simulations of hard spheres and disks are used to estimate the most probable cavity 
size, ^cavity) and a "rattle" size, prattle, over which a particle can translate holding all other par- 
ticles fixed. Both of these measures of free volume appear to extrapolate to zero at the random 
close packed density, p^cp, close to the density where extrapolations of the viscosity diverge. We 
also identify the onset of caging as the density at which cavities cluster. These results suggest 
that vitrification in hard core fluids can be viewed as a geometrical phenomenon, and that local 
free volume measures can identify the location of the onset of liquid-like dynamics, the complex 
dynamics of caged liquids, and vitrification. 
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Ensembles of hard spheres have long been utilized as prototypical models for liquids 
and solids. Since these systems can be defined purely by geometry, it seems reasonable 
to understand them in terms of void or "free" volume. Thus, a class of models, termed 
"free volume" models have been very popular in describing the transition of hard core fluids 
from the gas to the liquid state, from the liquid to the crystal and from the liquid to the 
glassy state 

mm 

. Inspite of the popularity of these ideas, the molecular connection 
between "free volume" and dynamic transitions in these fluids remains unclear ^. Previous 
molecular dynamics simulations have identified many dynamic signatures of the 

vitrification of hard core fluids and it is our goal to enunciate any underlying structural 
indicators. 

We simulate hard spheres and hard disks using the Monte Carlo method and periodic 
boundary conditions. The lateral size of the system, L, and the number of particles, A^, 
are held fixed in each simulation. We used L=100 in two dimensions {D=2) and up to 
L=60 for D=3. Density, or equivalently volume fraction, is defined as p = in D=2 

and p = ^^"p' D=3 where a is the particle diameter. (We report lengths in units of 
cr.) We generate random initial configurations, and thus some particle overlap was initally 
unavoidable at the highest densities. Equilibration was monitored by comparing various 
system properties, including the pressure and the chemical potential, to literature estimates 
(for example Q]). We also confirmed that system did not "age" as a function of time over 
the course of the simulations (10^ attempted moves/particle at all densities). 

Representative particle configurations for relatively low (p=0.40) and high density 
(p=0.53) hard disc {D=2) fluids are shown in Fig. ^ Note that although the disks are 
themselves of diameter 1, we have drawn regions of exclusion of diameter 2 in this figure. 
Thus, the unshaded areas correspond to regions where centers of new particles can be ac- 
comodated in the fluid without overlap. We observe that the "voids" in this figure become 
smaller and sparser as the p increases, and it is this effect we quantify in this paper. 

We evaluate two characteristic structural size scales using algorithms based on particle 
insertions . Following Ref . , "cavity" size distributions are defined by randomly 

selecting points in space and determining the size of the largest spherical cavity that can be 
inserted about them. This definition of a "cavity" is notably different from that employed 



in refs. 



14| . The distribution function for these maximal spheres displays a single 
peak, which is the most probable cavity size, ^cavity We expect ^cavity to be related to 



2 



characteristic cavity sizes measured by positron annihilation hfetime spectroscopy, but the 
precise relationship has not been conclusively demonstrated. We find that ^cavity varies 
linearly with Since both and the viscosity of these fiuids are expected to diverge 

near the random closed pack limit Q, Q, Q, Q|, a phenomenological connection 
between system dynamics and the presence of cavities is implied. We note that, although 
the most probable cavity size is expected to go to zero in this limit, cavities are still expected 
to exist beyond this point. 

Another useful measure of free volume is the "rattle volume" which corresponds to the 
volume explored by a particle given the constraints of surrounding particles [J, Q, ^ . Rat- 
tle size is evaluated in two equivalent ways. We consider the equilibrium distribution of 
acceptable displacements, Pdisp{r), made by randomly selected particles while holding all 
other particles fixed (see Fig. ^p). The maximum test displacement size is the box size so 
that all possible displacements are sampled. Alternatively, we can obtain the same results 
by considering successful insertions of pairs of spheres of size unity into the fiuid without 
overlap, and constructing a probability distribution function for their separations. Notably, 
Pdispi'f') has a peak at small r values for p ^ 0.25 in D=3, and for p ^ 0.35 in D=2 (see 
Fig. inset). Since the chosen particle is hindered by surrounding particles for distances 
)eyond the peak, we identify the peak as the most probable rattle size, prattle- Past work 
^ suggests that this rattle size is closely related to the Debye- Waller factor derived from 
dynamic scattering experiments. 

The fact that there is no peak in pdisp(^) for p < 0.25 (-0=3) and p < 0.35 {D=2), 
implies that there is no well-defined rattle volume in this density range, i.e., particles can 
freely move over the whole volume, as expected for a gas. For higher densities, however, our 
results suggest that rattle sizes are no longer "extensive" (i.e., scales with the system size), 
suggesting that holes are localized [see Fig. [T]. This is consistent with the fact that p ~ 0.26 
has previously been assigned as the onset of the "fiuid" region for hard sphere dynamics 
{D=3), where the free volume no longer percolates j^. Hence, prattle identifies the region 
where hard fiuid transport becomes liquid-like. Note that this density is much less than the 
density at the onset of "caging" in the dynamics, i.e., p = 0.53 for D=3, as we shall discuss 
below Q]. To further test the utility of prattle, "we consider the Lindemann criterion as 
applied to an amorphous phase [2^. This criterion requires that the ratio of the mean rattle 
distance to the interparticle distance be ~ 0.1-0.15 at the transition density. Specifically, for 
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both spheres and disks we find that .^rattle P^'''^ ~ 0.08±0.02 at the crystalhzation density, 
consistent with this generahzed Lindemann criterion. These two resuhs support the utihty 
of these free volume ideas. 

Fig. 121 shows our results for .^rattle for a range of p in both D=3 and D=2. We also 
plot the experimental cage sizes obtained from the dynamics of hard spheres to demonstrate 
the connection between equilibrum properties and dynamics in these situations 0|. 
As a first step towards understanding the density dependence of prattle, we note that in a 
simple "cell" model of a dense homogeneous fluid, the cage size, ^, is comparable to the 
average interparticle distance, so that ^ scales as p^^^^ at low densities and goes to in 
the limit of close packing. In Fig. El we compare prattle to the interpolation formula, ^ oc 
^p^^ / _ connecting these two limits. We employed this formula with proportionality 
constants a;=0.87 (D=2), and a=lA (-0=3) and literature values of prcp= 0.83 (-0=2) and 
0.644 (D=3), respectively. These results are consistent with the notion that this molecular 
measure of free volume vanishes in the limit of random closed packing. Alternate forms 
often employed for studying dynamic properties, such as prattle ~ (Prcp/p — 1)'^ or prattle ~ 
(1 — p/ PrcpV 1 also provide adequate fits to these data, reiterating this fact. The cavity size, 
^cavity, also follows the same density dependence with a=0.345 (-0=3) [Fig. E]. However, for 
-0=2 we had to fit both a=0.038 and prcp=0.79 to obtain good agreement with simulations, 
due to the large error bars associated with this small quantity. Note that the number of 
data points for ^cavity are typically different from prattle due to sampling issues that make the 
evaluation of prattle difficult at high densities. 

These plots exhibit the usual problem that the value of p where prattle and ^cavity 
depends sensitively on how we extrapolate to higher densities. To illustrate this point, we 
note that both .^rattle and .^cavity vary nearly linearly with p for 0.3< p < 0.45 for .0=3, 
and for 0.5< p < 0.65 for -0=2. Curiously, linear extrapolations of the low density prattle 



and .^cavity data to zero yield pg=0.56±0.03 for .0=3 and pg=0.74±0.05 for -0=2 



221. The 



density estimate in i5=3 closely matches the reported experimental glass transition density. 
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Pg ~ 0.56— 0.58 presented by jS, 



, where diffusion coefficient data were extrapolated to 
zero over the restricted range p < 0.52. The significance of this coincidence is unclear, since 
experiments and simulations where higher density data are considered (e.g., jl5l.ll6. 17 . isj ]) 
suggest that the viscosity does not diverge at pc, but rather at a density close to prcp- 
Similarly, we find that a plot of the logarithm of literature diffusivity data [l6| varies linearly 
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with ^r^ttie' agreement with the notion that random closed packing is close to the point 
where the system "vitrifies" jjj. However, we note that the determination of the exact 
density at which the viscosity diverges is somewhat uncertain due to the inherent errors in 
the high density data and the extrapolation process. For two dimensional disks, the "glass 
transition" density is harder to locate since most simulations introduce polydispersity to 
prevent crystallization. For polydisperse disks there is the estimate pc ~ 0.8 |2J, |25|, but 
this number is probably larger than for the monodisperse limit. Our pc estimate thus does 
match the experimental estimates obtained from extrapolation of low density data. 

Above we noted that visual examination of fluid configurations suggest that small, com- 
plex shaped cavities become increasingly separated with increasing p (Fig. We next 
define a scale relevant to the typical separation of the "void clusters" . We expect that spa- 
tial correlations associated with these void structures may have some relation to dynamic 



heterogeneities in the dense fiuid regime |2J, |2a, |22l , but the exact relationship is currently 
unclear. Specifically, we consider the distance away from a given insertion point that we 
must go such that the probability of insertions in a spherical (circular in D=2) shell is larger 
than the peak probability of insertion, located in the immediate vicinity (r < 1) of the first 
insertion point. For example. Fig. shows that this "crossover" scale, Across, at which 
insertions become more probable than in the immediate neighborhood occurs at ?5i 8 for 
p = 0.4. For randomly chosen insertions points, we expect the scale for successive insertions 
to be ^rand = {^/PinsY^^ , which uaturally increases for increased p. [pms is the insertion 
probability for a sphere into a snapshot of the fiuid.] It is natural to then consider the role 
of correlations by comparing the behavior of Across and ,^ran- Fig. HI shows both Across and ^^an 
as a function of p. At small p, .^cmss < ■Cran, apparently refiecting the tendency of particles 
to be clustered in the fiuid (Fig. Q). There is a characteristic p where Across first exceeds 
,^ran, and under these conditions we suggest that it is the holes that now cluster. Notably, 
the density where Across = ^ran corresponds to the density identified as the onset of "caged" 



dynamics in both D=2 (p ^0.48) 28] and D=3 (p ^0.53) |19|,|23|. Hence, we tentatively 
identify the point where the holes cluster as the onset of complex, structured fiuid behavior, 
and possibly a structural indicator of the onset of dynamical heterogeniety. 

It is useful to contrast our work to Ref. j^], which has shown that the rattle volume van- 
ishes in the vicinity of the Vogel temperature for Lennard- Jones particles. For hard spheres, 
our results suggest that random close packing plays the role of the "Vogel" point, i.e., the 
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density where the viscosity diverges and the free volume extrapolates to zero. Additionally, 
an alternate measure of the range of structure provides an estimate of the onset of caging. 

We acknowledge the National Science Foundation , Division of Materials Research for 
funding, and Ralph Colby, Pablo Debenedetti and Robin Speedy for useful discussions. 
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Figure Captions 



1. (a) Typical snapshots of the D=2 fluid. We siiow small portions of the simulation cell 
for ease of visualization, and the circles drawn have a radius of 1. Thus, the white 
spaces in these flgures represent areas where centers of additional disks can be placed 
without overlap. The left snapshot is for a density p=0.4, while the right picture is at 
p=0.53. (b) Main plot: Pdisp('")/47r as a function of r for the three dimensional hard 
sphere fluid at p=0.4. The two arrows represent characteristic sizes: the one closest 
to the origin is prattle, while the other distance is Across- The Pdispij) is used to generate 
y{r) = ^^^^wjT^S=T^ 1 where r(a;) is the gamma function, and 2'K^/'^r^~^ /T{r)dr is the 
volume of a D-dimensional hypersphere shell of radius r and thickness dr. The distance 
and density dependence of y{r) for p <0.495 for D=3 are consistent with past results 
j^. Inset: Plot of Pdisp('")/4vr as a function of r for the three dimensional hard sphere 
fluid at a variety of densities as sketched in the flgure. 

2. Plots of prattle a fuuctiou of p for D=2 (squares) and D=3 (circles). Uncertainties 
correspond to standard deviations in this data as obtained from block averages. The 



triangles are experimental cage sizes for hard spheres as reported by j21[ . Here we have 
normalized the experimental cage sizes by the bead diameter. Since our deflnition of 
prattle niay bc expected to be smaller than the cage size by a factor of 2, we divide the 



experimental data by 2 to ensure a proper comparison. The full line is the best flt to 
the cell model form prattle = « (^^) ^ — 1 
lower density data as discussed in the text. 



while the dotted line is a linear flt to 



3. Cavity size, ^cavity plotted as a function of density. Main plot is for hard spheres and 
the inset is for disks. Uncertainties are standard deviations in the data. The dotted 
lines are best flts over the low density ranges as discussed in the text, while the full 



line is the best flt to the cell model form ^cavity = ot 



4. Data for the Across (squares) and ^ran (lines) for hard spheres [main plot] and disks 
[inset] as a function of density. 
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